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Abstract. The 3:2 spin-orbit resonance between the rotational and orbital motions of Mercury results from a 
functional dependance on a tidal friction adding to a non-zero eccentricity with a permanent asymmetry in the 
equatorial plane of the planet. The upcoming space missions, MESSENGER and BepiColombo with onboard 
^ ' instrumentation capable of measuring the Mercury's rotational parameters, stimulate the objective to attempt to 

oo : an accurate theory of the planet's rotation. We have used our BJV relativistic model of solar system integration 

including the spin-orbit motion of the Moon. This model had been previously built in accordance with the 
requirements of the Lunar Laser Ranging observational accuracy. We extended this model to the spin-orbit 
couplings of the terrestrial planets including Mercury; the updated model is called SONYR (acronym of Spin-Orbit 
N-BodY Relativistic model). An accurate rotation of Mercury has been then obtained. Moreover, the conception of 
, the SONYR model is suitable for analyzing the different families of hermean rotational librations. We accurately 

identify the non- linear relations between the rotation of Mercury and its dynamical figure (C/MR^, C20, and 
C22). Notably, for a variation of 1% on the C/MR^ value, signatures in the ip hermean libration in longitude as 
well as in the rj obliquity of the planet are respectively 0.45 arcseconds (as) and 2.4 milliarcseconds (mas). These 
^ ' determinations provide new constraints on the internal structure of Mercury to be discussed with the expected 

^ I accuracy forecasted in the BepiColombo mission (respectively 3.2 and 3.7 as according to Milani et al. 2001). 
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^ 1. Introduction sertions introduced by Peale in 1976. According to this 

' author, the determination of the four parameters, namely 

. - . Pettengill and Dyce discovered in 1965 the 3:2 spin-orbit ^^^^ ^^^^ ^ ^^^^^^ sufficient to constrain size 

resonance state of Mercury (the rotational and orbital pe- ^^^^^ ^j^^ Mercury's core (C20 and C22 are spher- 

riods are 56.646 and 87.969 days respectively). This partic- .^^^ harmonics of the second degree, is the obliquity of 

ular resonance results from a functional dependance on a Mercury, and ^ is the libration in longitude of 88 days), 
tidal friction adding to a non-zero eccentricity with a per- 
manent asymmetry in the equatorial plane of the planet. 

Besides, in a tidally evolved system, the spin pole is ex- In the first section of this paper, we present our ap- 
pected to be trapped in a Cassini state: the orbital and proach which takes into account the coupled spin-orbit 
rotational parameters are indeed matched in such a way motion such as it derives from the integration of the en- 
that the spin pole, the orbit pole, and the solar system tire solar system. In the second section, we present the 
invariable pole remain coplanar while the spin and orbital dynamical behavior of the Mercury's spin-orbit motion, 
poles on average precess at the same rate (Colombo 1965; In the last section, following our previous analysis of the 
Peale 1969, 1973). different families of hermean rotational librations (see our 
The upcoming missions, MESSENGER (Solomon previous paper Rambaux and Bois 2003), we make in evi- 
et al. 2001) and BepiColombo (Milani et al. 2001) with ^ence some accurate signatures due to parameters related 
onboard instrumentation capable of measuring the rota- to the dynamical figure of Mercury (C/MR^, C20, and 
tional parameters of Mercury stimulate the objective to C'22) in the librations as well as in the obliquity of the 
attempt to an accurate theory of the Mercury's rotation, planet. We find 0.45 arcseconds (as) and 2.4 miUiarcsec- 
The current method for obtaining constraints on the state o^ds (mas) respectively on ip and 77 for a variation of 
and structure of the hermean nucleus is based on the as- 1% on the C/MR^ value. These determinations provide 
new constraints on the internal structure of Mercury to 
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BepiColombo mission (respectively 3.2 and 3.7 as accord- 
ing to Milani et al. 2001). 

2. The Spin-Orbit N-Body Relativistic model 

For obtaining the real motion of Mercury, we have used 
our BJV model of solar system integration including the 
coupled spin-orbit motion of the Moon. This model, ex- 
panded in a relativistic framework, had been previously 
built in accordance with the requirements of the Lunar 
Laser Ranging observational accuracy (Bois 2000; Bois & 
Vokrouhlicky 1995). We extended the BJV model by gen- 
eralizing the spin-orbit couplings to the terrestrial plan- 
ets (Mercury, Venus, Earth, and Mars). The model is at 
present called SONYR (acronym of Spin-Orbit N-BodY 
Relativistic model). As a consequence, the SONYR model 
gives an accurate simultaneous integration of the spin- 
orbit motion of Mercury. The integration of the solar sys- 
tem, including the Mercury's spin-orbit motion, uses a 
global reference system given by the solar system barycen- 
ter. The model is solved by modular numerical integra- 
tion and controlled in function of the different physical 
contributions and parameters taken into account. It per- 
mits the analysis of the different families of librations and 
the identification of their causes (such as the planetary 
interactions) or their interdependences with the param- 
eters involved in the dynamical figure of the planet. A 
complete description of the SONYR model is given in our 
previous paper (Rambaux and Bois 2003). The spin-orbit 
motion of Mercury is characterized by two proper frequen- 
cies (namely $ = 15.847 years and = 1066 years) and 
its 3:2 resonance presents a second synchronism, which 
can be understood as a spin-orbit secular resonance (H = 
278 898 years). 

Using the SONYR model and its method of analysis 
(modular and controlled numerical integration, differenti- 
ation method), the librations are accurately surrounded 
and identified. Our previous paper Rambaux and Bois 
(2003) presents a detailed analysis of the whole spin-orbit 
mechanism of Mercury as well as its main librations. 

3. The Mercury's spin-orbit motion 

Using our SONYR model, the present work is devoted for 
studiing the impact of the dynamical figure of Mercury on 
its spin-orbit motion. We have then analyzed the impact 
of the three dynamical parameters C/MR^, C20, and C22 
on the rotational motion of Mercury. 

Figure Q presents the dynamical behavior of the 
Mercury's rotation expressed by the Euler angles -0, 6*, and 
ifi (three first panels) as well as its obliquity 77 plotted in 
the fourth panel over 500 days. The Euler angles (^/j, 0, 1^) 
related to the 3-1-3 angular sequence describe the evo- 
lution of the body-fixed axes Oxyz with respect to the 
axes of the local reference frame OXYZ . Let us recall the 
definition used for these angles : ip is the precession an- 
gle of the polar axis Oz around the reference axis OZ, 
is the nutation angle representing the inclination of Oz 
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Fig. 1. The rotational motion of Mercury expressed in the 
classical 3-1-3 Eulerian sequence ('0, 9, f) with respect to 
the ecliptic reference frame. In the bottom panel is plotted 
the obliquity 77. The four curves are plotted over 500 days 
(horizontal axes). Degrees are on the vertical axes for ■!/;, 0, 
and iy9, while the vertical axis for r] is plotted in arcminutes 
(amin). 



with respect to OZ, and (p is the rotation around Oz and 
conventionally understood as the rotation of the greatest 
energy (it is generally called the proper rotation) . The axis 
of inertia around which is applied the proper rotation is 
called the axis of figure and defines the North pole of the 
rotation (Bois 1992). Let us remark that in this Figure, 
we have removed the mean rotation of 58.646 days in the 
Lp angle in order to better distinguish the librations. The 
resulting amplitude of libration appearing in Lp is 20 as, 
which is in good agreement with Balogh and Giamperi 
(2002). 

The integration of the SONYR model including the 
simultaneous spin-orbit motion of Mercury gives also the 
dynamical behavior of the hermean obliquity by the way 
of the following relation : 



cos rj = cos i cos 9 + sin i sin 6 cos(f2 — 0) 



(1) 



where i and are respectively the inclination and ascend- 
ing node of the Mercury's orbital plane relative to the 
ecliptic plane (the obliquity is the angle between the po- 
lar axis and the orbital plane). In our computations pre- 
sented in this paper, the used initial conditions are listed 
in Tabled For 77, 0, and 0, we use the values determined 
in our previous work (Rambaux and Bois 2003). 

The dynamical structure of the Mercury's rotation 
given in Figure^results from the direct solar torque whose 
coordinates include the planetary interactions. The effect 
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Fig. 2. Signatures of the variations (with a step equal to 0.1%) of the dynamical parameter values on the Mercury's 
obliquity: (a) C/MR^, (b) C20, and (c) C22. On each panel, 10 curves are plotted over 500 days; the maximum of 
amplitude is obtained for 1%. Milliarcseconds (mas) are on the vertical axes. 



of the gravitational torque due to Venus is five orders of 
magnitude smaller than the one due to the Sun. The 
rotation period of 58.646 days appears in the ip and 9 an- 
gles, as well as in 77. Signature of the Pa orbital period of 
87.969 days is clearly visible in the ip angle (this angle is 
called in literature the libration in longitude of 88 days). 
A third period appears in the ipj and 77 angles, namely 
176.1 days. This period P results from the 3:2 spin-orbit 
resonance {P = 2P\ = 3P^). 

Table 1. Our initial conditions at 07.01.1969 (equinox 
J2000): (a) Mean values derived from the SONYR model; 
(b) Seidelmann et al. (2002). In addition, our SONYR 
value of the Mercury's obliquity is 1.6 amin. 

Mercury 
Rotation angles 



e 



48.386 deg 1/; 
7.031 deg 
299.070 deg ip 



= -0.616 lO"'^ deg/day (a) 
e = -0.267 10"** deg/day (a) 
= 6.138505 deg/day (b) 



4. The impact of the dynamical parameters 

The gravity field of Mercury is globally unknown. The 
tracking data from the three fly-by of Mariner 10 in 1974- 
75 have been re-analyzed by Anderson et al. (1987) in 



order to give an accurate estimation of the normalized 
coefficients C20 and C22. The resulting nominal values are 
C/MR2=0.34, C20 = (6.0 ± 2.0) 10-^ and C22 = (1.0 ± 
0.5) 10~^. As a consequence, the dynamical coefficients 
a = (C - B)/A, (3 = [C - A)/B, and 7 (B - A)/C 
(where A, B, and C are the three principal moments of 
inertia of Mercury) are infered from C20 and C22 by the 
following formulae : 



-Co. 



2C2 



7 



C20 - 2C22 + C/MR? 
—C20 + 2C22 



C20 + 2C22 
4C22 



C/MR? 



(2) 



In these conditions, we may assume the uncertainties on 
the Mercury's dynamical parameter values as dynamical 
variations on the Mercury's rotation. We have computed 
the impact of the variations of the greatest principal mo- 
ment of inertia, C/MR^, as well as C20 and C22, both on 
the obliquity and the libration in longitude. Figures Eland 
13 present the results plotted over 500 days. 

These Figures show three panels corresponding to sig- 
natures of the variations of C/MR^ (a), C20 (b), and C22 
(c) on 7] and ip. In each panel, 10 integrations with a step of 
0.1% over the given parameter are plotted over 500 days; 
the maximum of amplitude is obtained for 1%. In each 
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Fig. 3. Signatures of the variations (with a step equal to 0.1%) of the dynamical parameter values on the Mercury's 
libration in longitude: (a) C/MR^, (b) C20, and (c) C22- On each panel, 10 curves are plotted over 500 days; the 
maximum of amplitude is obtained for 1%. Arcseconds (as) are on the vertical axes of the top and bottom panels 
while it is about milliarcseconds (mas) in the middle panel. 



plot, the maximum amplitude is reached with the charac- 
teristic period of 176.1 days for the obliquity 77, and 87.969 
days for the libration in longitude ip. The results related 
to variations of 1% are listed in Table El 

Table 2. Results of the 1% variations of the dynamical 
parameters on the Mercury's obliquity and libration in 
longitude. 





C/MR-' 


C20 


C22 


V 


2.4 mas 


1.8 mas 


0.7 mas 


f 


0.45 as 


14 mas 


0.45 as 



In the Figure 12 a variation of 1 % on the C /MR^ pa- 
rameter induces a peak to peak amplitude of 2.4 mas in the 
obliquity, whereas with the C20 and C22 coefficients, the 
amplitudes are 1.8 and 0.7 mas respectively. The charac- 
teristic period of these signatures is P = 176.1 days, which 
expresses the spin-orbit coupling of Mercury. The dynami- 
cal behavior of 77 proves to be the right tracer for analyzing 
the relationship between the rotation of Mercury and its 
dynamical figure. 

In the Figure O the impact of variations of the C20 
value on the libration in longitude is 30 times smaller than 
the one due to variations of C /MR^ and C22 . This is make 



clear by the fact that the dynamical parameter 7 depends 
only on C22 and C/MR^. As a consequence, C20 has only 
an indirect action on the if angle. In addition, the dynam- 
ical behaviors in FigureOl panels (a) and (c), are symmet- 
rical; this fact reflects structure of the third equation in 
the formulae (2). The relation between C20 and proves 
to be a good tracer of the P spin-orbit period while the 
relations between C/MR^ and or C22 and (/s, express 
the Pa period. 

5. Conclusion 

One important objective of the upcoming space missions, 
MESSENGER and BepiColombo, is to constrain the na- 
ture of the Mercury's core. The theoretical framework fol- 
lows the procedure proposed by Peale (1976) and Peale 
et al. (2002). The onboard instrumentations are planed 
for measuring the Mercury's rotational parameters with 
a great accuracy. However, our work shows that for ob- 
taining an excepted accuracy up to 1% on the C/MR^ 
coefficient, the measurements have to reach at least 0.45 
as on the Mercury's rotation and 2.4 mas on the obliq- 
uity. To our knowledge, these high accuracies are not yet 
achieved within the simulations. 

The Mercury's rotation strongly depends on the 
C/MR2 (0.45 as) and C22 (0.45 as) coefiicients. Besides, 
the obliquity of Mercury hierarchically depends on C /MR^ 
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(2.4 mas), C20 (1-8 mas), and C22 (0.7 mas). The dynami- 
cal behavior of Mercury's obliquity is finally the tracer for 
analyzing the non-linear relations between the rotation of 
Mercury and parameters of its dynamical figure. 
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